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Abstract 

This paper presents a novel vibration-based piezoelectric energy harvester capable of passively 
tuning its resonant frequency to a wide range of frequencies. The device comprises a dual 
bimorph with a mass at its free end. A novel sliding mechanism, consisting of two oblique 
springs connected to the tip mass, is proposed to widen the resonance frequency of the device 
even to very low frequencies. The application of two oblique springs results in an additional 
stiffness and axial load that are introduced within the system, such that the resonance 
frequency of the device is now a function of both the stiffness and axial load associated with 
the spring forces. An operator can manually change the resonance frequency of the harvester 
just by small adjustments of the sliding mechanism. Further, the device allows one to tune the 
resonance frequency of the beam to match very low frequencies without the requirement of 
having a large proof mass. The analytical solution of an axially loaded cantilevered 
piezoelectric energy harvester with tip stiffness, using Euler-Bernoulli beam assumptions, is 
presented. A parametric case study is presented to demonstrate the performance of the device. 

(Some figures may appear in colour only in the online journal) 



1. Introduction forms of piezoelectric energy harvesters. The main drawback 

of these devices is that to be efficient, their resonance 

Vibration-based energy harvesting has become the focus frequency needs to match exactly the frequency of the ambient 

of considerable research effort in the last decade. This vibration [2-8]. Even small frequency mismatches between 

increase in research has been brought on by modern the harvester resonant frequency and the dominant frequency 

advances in low-power-consumption devices and the need of excitation can cause significant decreases in power output. 
to provide wireless technology. Vibration energy can be Several approaches have been proposed to broaden the 

converted to electrical energy through three main kinds frequency range of these devices, such as devices with 

of transduction mechanism: electrostatic, electromagnetic multiple cantilever beams with different lengths and tip 

and piezoelectric. Among these transduction mechanisms, masses [3]. It has been shown that multiple beams can 

piezoelectric transduction has recently received the greatest efficiently operate in a wide band of frequencies; however, 

attention mainly due to the large power density, ease of the geometry of this kind of device may negatively affect 

application and direct conversion of mechanical to electrical the cost and size of the harvester. A nonlinear mechanical 

energy [1]. resonator has been proposed in several studies to enlarge 

Cantilevered beams with one (i.e., unimorph) or two the energy harvester's power bandwidth [9-13]. In Cornwell 

(i.e., bimorph) piezoceramic layers are the most popular et al [14] a tuned axillary structure was attached to the host 

structure to overcome the limitation of low power output at 

1 Author to whom any correspondence should be addressed. a low excitation frequency; however, this method requires 
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Figure 1. The proposed tunable energy harvester: (a) 3D view and (b) plan view . 



a prior knowledge about the excitation frequency of the 
host structure. An active self-tuning device was proposed 
in [15]. The main limitation of an active approach is the 
requirement of additional power to change the dynamic or 
electronic properties of the harvester [16]. In contrast, a 
passive approach only requires an initial power to tune the 
device; as a result, interest in this topic has continually 
grown in the past few years. Wu et al [17] proposed a 
movable mass to change the frequency of the resonator. 
A self- tunable piezoelectric harvester of movable cantilever 
couples was investigated by Jo et al [18]. The proposed 
device was able to work in two phases in which four tuned 
frequencies could be achieved. Gu and Livermore [19, 20] 
exploited the centrifugal force of a rotational system to 
provide tensile stress and change the resonance frequency 
of a harvester for rotating applications. A novel simply 
supported bimorph with a variable compressive axial preload 
was proposed in [2] in order to reduce the stiffness of the 
beam and, as a consequence, the resonance frequency of the 
harvester. It has been shown that the compressive load can be 
used to tune the harvester and to improve the efficiency of 
the piezoceramic layers by increasing the electromechanical 
coupling coefficient [21]. An axial preload effect on the power 
generation and shift frequency of a piezoelectric cantilevered 
beam was numerically investigated by Hu et al [22]. Eichhorn 
et al [23] proposed a piezoelectric cantilevered beam with 
two additional arms that were used to provide compressive 
and tensile load to the tip of the beam in order to reduce and 
increase the frequency of the energy harvester. A magnetic 
force was used in [4] in order to change the tip stiffness and as 
a result the frequency tuning of a cantilever beam; however, 
the performance of this device in an environment with 
magnetic interference may be not reliable. A low-frequency 
zigzag piezoelectric energy harvester was proposed in [24] 
and validated by a comprehensive parametric study. 

In this work a novel tunable vibration energy scavenger is 
proposed that provides mechanically variable and adjustable 
tip stiffness and axial load for a bimorph piezoelectric 
energy harvester. An analytical model is used to show the 
performance of the device. 




Figure 2. The schematic geometry of the device. 



The paper begins by an overview of the proposed device. 
An analytical solution for the vibrating bimorph cantilever 
beam with tip stiffness and axial load is presented in section 3. 
This is followed by section 4, which reports the results of a 
case study, followed by the conclusion in section 5. 

2. Overview of the proposed tunable energy 
harvester 

The proposed device consists of two parallel bimorphs 
clamped at one end and connected to each other at the other 
end, as shown in figure 1. A sliding mechanism is used to 
change the equivalent vertical stiffness and the axial load 
applied by two springs attached to the tip of the dual beam. 
As a result, the dynamic behavior of the system can be tuned 
manually by an operator according to the dominant frequency 
of the intended environment. 

A kinematic mechanism associated to a given displace- 
ment (d) of the sliding part is shown in figure 2. Theoretically, 
two extreme positions can be defined, that is, a = 0° and 
90°, which correspond to minimum and maximum stiffness, 
respectively. However, it is recommended to impose an initial 
position (Go or ofo) for the mechanism so that it can be changed 
around this, to increase or decrease the resonance frequency 
of the energy harvester, as will be shown in section 4. Based 
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N=2 N s cos(0) 
K=2K s sin 2 (9) 



Figure 3. A schematic view of the device with its equivalent cantilever beam with mass, axial load and stiffness at the tip. 



a) 

M t : Tip mass 

K: Tip stiffness 

TV: Tip axial load 

Rf. Circuit load 

h p : Piezoelectric thickness 

h' Substructure thickness 
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Figure 4. Configurations of a geometrically uniform bimorph piezoelectric energy harvester with tip axial load and stiffness and excited by 
translation and rotation of its base, (a) Circuit connection of the piezoceramic layers of bimorph with relative displacement and (b) their 
cross- sectional view. 



on geometrical considerations and since L\ and e are fixed 



a — arccos 



6 — arctan 



ao±d 

u 

L\ sin (a) 



(i) 

(2) 



where ao, d, e and L\ are shown in figure 2. 

To investigate the dynamic behavior of the proposed 
device the equivalent beam shown in figure 3 is considered. 
Under the assumption of small displacements, the vertical 
stiffness (K) and axial load (N) at the tip of the beam 
are proportional to sin 2 (£) and cos(0), respectively [25]; 
therefore, the following relationships can be defined: 



K = 2K s sm 2 (0) 
N = 2N s cos(0) 



(3a) 
(3b) 



where ^ s is the stiffness of each spring, is the spring angle 
to the x axis of the beam and N s is the force exerted by each 
spring. It can be observed that the values of K and N can be 
modified by changing the angle 9. The damping component 
as a result of the tip spring is not shown in figure 3; however, 
in section 3 it will be considered in the equation of motion. 

In section 3 the equations governing the electromechani- 
cal behavior of the proposed device are presented. 



3. An analytical model of a cantilevered bimorph 
with tip mass, stiffness and axial load 

This section presents the modal analysis of the cantilevered 
bimorph shown in figure 4. To the best of the authors' 
knowledge, the dynamic response analysis of a piezoelectric 
cantilevered beam with tip stiffness and axial load has not 
been reported in any previous research. 

By having the modal frequencies and shape functions in 
the same fashion as presented in [26], the voltage and power 
frequency response function of the system will be defined in a 
closed-form solution. 

3.1. Coupled mechanical equation and modal analysis 

The cantilever beam shown in figure 4 is subjected to arbitrary 
translation and rotation motions of its base, which are denoted 
by g(t) and h(t) respectively. The effective transverse base 
displacement Wb(x, t) can be written as [26] 



Ub{x,t) =g(t)+xh(t). 



(4) 



Using the Euler-Bernoulli beam theory, the partial differential 
equation for an axially loaded cantilever beam with tip mass 
under base excitation has been defined by Erturk et al [26]. On 
the addition of tip axial load the partial differential equation 
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is [27] 



d 2 M{x, i) 



+ cJ 



9 WrelO, t) 



dx 2 

+ c t 8(x — L c ) 



dx 4 dt 

du K \(x,t) 
dt 



+ m 



d U re \(x, t) 3 U xt \{x, t) 



dt 2 



+ N- 



= -[m + M t 8(x-L c )] 



dx 2 
9 2 Wb(x, t) 

dt 2 



(5) 



where M(x, t) is the internal bending moment, c s is the 
equivalent damping coefficient of the composite cross section, 
c t is the equivalent vertical damping coefficient of the springs 
at the tip of the beam, / is the equivalent area moment of 
inertia of the composite cross section, m is the mass per 
unit length of the beam, L c is the length of beam and 8 (x) 
is the Dirac delta function. It should be noted that N is not 
tangential, that is, the direction of N is considered to be always 
horizontal. 

In equation (5) m is defined as 



m = b(p s h s + 2/>p/z p ) 



(6) 



where b is the width, p s is the mass density of the substructure 
material, p p is the density of the piezoceramic material, 
h s is the substructure thickness and h p is the piezoelectric 
thickness, as shown in figure 4(b). 

Instead of defining damping coefficients in the physical 
equation of motion, an undamped system is considered, by 
setting c s I and c t to zero in equation (5) [26]; therefore, 
the modal damping is used in the equation in modal 
coordinates, as is common in practice. It is noteworthy that 
it is common in many dynamic problems for the sake of 
simplicity to approximate the different damping components 
by an equivalent modal damping ratio. The internal bending 
moment M(x, i) is the integral of the first moment of the axial 
stress over the cross section [26] 



M(x, t) = b 



X 
f 



T?zdz + f 2 Thdz 

hs i / h s i 



T\zdz 



(7) 



where T\ and T\ are the axial stress in the piezoelectric and 
the main structure layer, respectively, which are given by the 
following constitutive relations: 



T\ 



y&, 



i l 



^E C P 



c n SK - e 3 iE 3 



(8) 



where Y s is the Young's modulus of the substructure layer, 
c^ is the Young's modulus of the piezoeceramic layer at 
constant electric field, e$i is the piezoelectric constant and 
E 3 is the electric field component, in the z direction. The 
subscript and superscript s and p stand for the substructure 
and piezoceramic layer respectively. Under the assumption 
of plane stress, the Young's modulus of the piezoeceramic 
layer can be expressed as cf. — 4- , where 5? x is the elastic 

compliance at constant electric field [28]. Also, e$i can be 



defined in terms of the piezoelectric constant d$i as e$i 



ds\ 



It should be noted that, since the electrodes covering the 
opposite transverse faces of the piezoelectric layers are very 
thin, as compared to the overall thickness of the harvester, 
their contribution to the thickness dimension is negligible. 
The axial strain components in the piezoelectric S\ and 
substructure S\ at location z from the neutral axis of the beam, 
based on Euler-Bernoulli beam theory, are proportional to the 
curvature as 



S\ (x, z, t) 



9 MrelOc, t) 

1 dx~ 2 * 



(9) 



The length to thickness aspect ratio of the beam is usually 
high enough to neglect the shear effect; as a result the 
Euler-Bernoulli theory can be considered an appropriate 
assumption. Equations (8) and (9) are obtained under the 
assumption that the layers are perfectly bonded and exhibit 
linear elastic material behavior. 

Considering that the piezoelectric layers of a bimorph can 
be combined in two configurations, that is, series or parallel, 
the electric field E 3 in equation (8) should expressed in terms 
of the respective voltage term (i.e., v s (0, v p (0) for each 
configuration [29] using the following relationships: 



E 3 (f) = 



E 3 (t) 



v s (0 
2/zp 



for series 



for parallel. 



(10a) 



(10/7) 



The inertial bending moment can be obtained for a series or 
parallel circuit from equation (7) as [26] 



M(x, t) : 



-YI 



d 2 U re \(x,t) 

~dx~ 2 



+ Tv{t)[H{x)-H{x-L c )} (11) 



where YI is the bending stiffness term of the composite cross 
section for a constant electric field, T is the coefficient of 
the backward coupling term and H(x) is the Heaviside step 
function. 

Substituting equation (11) into (5) the coupled beam 
equation for the series configuration can be obtained as 

9Vel(*, 9 5 w re i(x, t) 

YI — V c s I- 



dx 4 
+ c t 5 (x — L c ) 



dx 4 dt 

dw re i(x, t) 
dt 



d 2 u vt \{x,i) d 2 u vt \{x,i) 

d8(x) d8(x-L c ) 



I> s (0 



dx 



[m + M t 8{x-L c )] 



dx 

9 2 Wb(x, t) 
dt 2 



(12) 



The coefficient of the backward coupling term for the series 
configuration in equation (12) is [26] 



e 3 \b 

2/Zn 



K + 



K 



h l 
4 



(13) 
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This coefficient for the parallel configuration is 2T. The 
coupled beam equation for parallel configurations can be 
obtained from equation (12) by changing F to 2I\ 

The equivalent bending stiffness term YI in equation (12) 
is [20] 



YI ■ 



2b 
~3 



hi 



+ tf 



K + 



h 



hi 



(14) 



To solve the differential equation (12) using the method of 
separation of variables [30], four boundary conditions are 
needed. 

The boundary conditions at x = are 



Krel(0, t) = 0, 



3w re i(x, i) 



dx 



0. 



(15) 



x=0 



The boundary condition at x — L c for the bending moment can 
be expressed as [20] 



YI 



d U re \(x,t) d 3 u re \(x, t) 



dx 2 



+ /t" 



dfldx 







(16) 



x=L c 



where 7 t is the mass moment of inertia of the proof mass; the 
shear boundary condition considering axial load is [31] 

9 3 w re i(x, t) d 2 u Y e\(x, t) 

YI \ - -Mi- 



dx 3 



KUrd(x,t)+N 



dt 2 

d 2 u K \(x,t) 

dx 2 



= 



(17) 



x=L c 



where K and N are defined in (3). Equation (17) is where the 
effect of tip stiffness shows itself. 

The method of separation of variables can be used to 
solve the differential equation defined in equation (12) as 



WrelfcO = <P(x)v(t) 



(18) 



where (p(x) is the modal eigenfunction and rj(t) is the 
modal mechanical coordinate expression. It should be 
mentioned that (p(x) is the same for both series and 
parallel configurations while r](t) is different. Based on 
the proportional damping assumption and substitution of 
equation (18) into equation (12) for free vibration leads to [27] 



YI d 4 <p(x) N d 2 <p(x) m d 2 r](t) 



= Y (19) 



4>(x) d 4 x <p(x) d 2 x 7}(t) d 2 t 

where y is a constant. Replacing y by mco 2 one can obtain two 
separated differential equations in the time and space domains 

as 



dVO . 2 M 







y/ d 4 0(x) | v d 2 0(x) 



■ mco (p(x) = 0. 



(20) 



(21) 



d 4 x d 2 x 

Equation (21) leads to frequency and mode-shape expressions 
for the free vibration, in which the axial load N is a basic 
dynamic parameter. Dividing equation (21) by YI leads to 



d 4 0(x) , 2 d 2 <p(x) 



d 4 x 



+ ^ 



d 2 x 



<2 4 0(X) = O 



(22) 



where q and a are defined as 



q = 



YI 



YI 



A general solution for equation (22) is [27] 



(23a) 
(23b) 



<p(x) — D\ sin(/?x) + D2 cos(^x) 

+ D3 sinh(sx) + D\ cosh(sx) (24) 

where s and fi can be obtained by solving the characteristic 
function of (s 4 + q 2 s 2 — a 4 = 0) for s, 



N 



q 4 



P = 



N 



4_l q 



(25a) 



(25b) 



Considering the boundary conditions in 
equations (15)— (17), the modal frequencies (eigenvalues) 
and the coefficients D\,D2,D^, and D\ can be obtained. 
Substitution of equation (24) into the geometric boundary 
conditions defined in equations (15)— (17), gives 



10 1" 


~Di 


P s 


D 2 


V31 A32 A33 A34 


O3 


Ul A42 A43 A44_ 


_D 4 



= 0. 



(26) 



In equation (26), the first two rows are derived from the 
boundary conditions defined in equation (15), whereas the Ay- 
are derived from the boundary conditions in equations (16) 
and (17) as follows: 



A31 



-kfio/ cos(fiL) - Yip 1 sin(fiL) 



A32 = hPo/ sm(fiL) 

A33 



YIp z cos(pL) 



hsco cosh (sL) + YIs sin h(sL) 



A34 = -I t sco 2 sinh (sL) + YI s 1 cos h(sL) 
A41 = M t co 2 sin (pL) - Yip 3 cos (pL) 

-K sin (pL) -Np cos (pL) 
A42 = M t co 2 cos(pL) + Yip 3 sin (pL) 

-Kcos(pL)+Npsin(pL) 
A43 = M t co 2 sin h (sL) + YIs 3 cos h (sL) 

— K sin h (sL) — Ns cos h (sL) 
A44 = M t co 2 cos h (sL) + YIs 3 sin h (sL) 

— K cos h (sL) — Ns sin h (sL) . 



(27) 
(28) 
(29) 
(30) 

(31) 

(32) 
(33) 
(34) 



For nontrivial solutions to exist, the following determinant is 
required to vanish, and this defines the characteristic function 
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to find modal frequencies or eigenvalues 

10 1 
P s 

A31 A32 A33 A34 
A41 A42 A43 A44 



20 



0. 



(35) 



Solving equation (35) for co, infinite eigenvalues co v are 
obtained. 

By substituting each eigenvalue co v in equations (23b) and 
(24), the mode shapes are obtained as 



<pi(x) = A r 



Pt 
sin (firx) sinh (e T x) 



where T r is 



+ T r (COS (fax) — COSh (8 t x)) 



k 



J; A 33r - A3i r 
A32r - A34 r 



(36) 



(37) 



and A r is a modal amplitude constant that can be evaluated 
using mass normalization based on orthogonality of the modal 
shape functions [29] 



/' 

Jo 



(p s (x)m(p r (x) dx + (p s (L c )M t (p r (L c ) 



~d0 s (x) d0 r (x) 
h 



dx 



dx 



= 5,, 



(38) 



where 8 sr is the Kronecker delta function. The normalized 
modal shape functions <j> T {x) are the same for series and 
parallel configurations based on the zero load (R\ = 0) 
assumption (short circuit). For non-zero values of the load 
resistance the modal coordinates Tj T (t) are different since a 
point moment as a result of voltage appears in the coupled 
beam equation (12) [26]. 

The beam has infinite vibration modes, so the general 
response is a linear combination of all vibration modes 



U Te \(x, t) = y^(p r (x)f] r (t). 



(39) 



It is noteworthy that one can find the critical load N cr by 
setting co = and solving equation (35) for N. 

To evaluate the behavior of this device, the dimensionless 
tip stiffness and axial load are defined as K = = and N = 



N 



respectively, where K 



3YI 



is the theoretical stiffness 



N ™ r .,, ~ ^ 

of a cantilever beam with the polynomial shape function 
assumption and TV = ^rM is the critical axial load for a 
cantilever beam without tip mass and stiffness. Figure 5 shows 
the stable and unstable regions for different combinations of 
K and N. It should be noted that as K approaches infinity N cr 
reaches the value of the critical axial load for a beam with 
clamped-pinned boundary conditions. 



<^ 



18 
16 
14 
12 
10|- 

8 

S 

4 

2 




Figure 5. The normalized critical axial compressive load for 
different tip normalized stiffnesses of a cantilevered beam. 



3.2. Electromechanical analysis of a bimorph with tip mass, 
stiffness and axial load 

A detailed discussion of the coupled electrical circuit 
equations of a piezoelectric layer for series and parallel 
configurations under dynamic bending is beyond the scope 
of this paper. Readers may refer to [26, 29, 32] for more 
details. Considering the piezoelectric constitutive equations 
and Gauss' law, the steady state voltage response of the 
bimorphs shown in figure 4 for a series configuration is [26] 



Vs(0 



E? 



)OJK x F v 



joj —(O+il'^COrOO 



R\ 



1 +j*4 + E? 



=1 ,£-„? 



}Q^rXr 



(40) 



+)2^(D X (D 



and for a parallel configuration is 



v P (0 = 



Er° 



)OJK x F x 



-(D+il'^xCDrQ) 



2Ri 



+j^c p + E" 



z l rn^—n 



Jft^rXr 



(41) 



'j0 2 +i2^ r (D r (D 

In equations (40) and (41), F T is the amplitude of the harmonic 
modal forcing function, co is the frequency of the base 
excitation, j is the unit imaginary number, £ r is the modal 
damping, co v is the modal frequency, /c r is the modal coupling 
term in the electric circuit, R\ is the resistive load in the electric 
circuit, C p is the capacitance of each piezoelectric layer and 
X r s and xF are modal electromechanical coupling terms. 

For the base harmonic excitation given by equation (4), 
in which g(t) = Yqq)^ 1 and h(t) = Woe?* *, the modal harmonic 
forcing function is F r e Jft ^, where F r is expressed as [26] 



F r = CO Z 



Jo 



mYo I r (x) dx + ra^o 



I X( 

Jo 



:0 r (x) dx 



+ M t r (L)(F o +M/ o ) 



(42) 



where Fo and ^0 are the translational and small rotational 
displacement amplitudes of the base. 
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The modal electromechanical coupling terms / r s and Xr 
for series and parallel are 

d0 r (x)~ 

x=L r 



Xi 



dx 



and 



xF = 



2r 



d0rW 

dx 



(43) 



x=L c 



where r is defined in equation (13). 

Also, the forward coupling term /c r is expressed as 



£31 



h pc b / 
JO 



d 2 r (x) 
d 2 X 



dx = e^ihpcb 



d0rW 
dx 



(44) 



x=L c 



where /z„ 



/*p+/*s 



As was mentioned before, as far as 



modal analysis is concerned, the bimorph is identical in 
series and parallel connection configurations. However, the 
backward piezoelectric coupling effect that is manifested in 
the tip moment induced by the coupling term k t for an 
arbitrary electrical load is different for series and parallel 
connections of the piezoceramic layers, which affects the 
vibration response of the cantilever (the last term on the left 
hand side of equation (12)) [33]. 

Under the assumption that the electrode pairs covering 
the top and bottom faces of the piezoceramic layers 
are perfectly conductive, a single potential difference can 
be defined across the electrodes. The internal (inherent) 
capacitances of the piezoelectric layers are considered in the 
circuit along with a resistive electrical load (denoted by R\ in 
figure 4). The piezoceramic layers are assumed to be identical 
and their conductive electrodes are assumed to fully cover the 
respective surfaces (top and bottom) of these layers. 

The capacitance C p for each layer is 



C P = 



-33 



bL 



K 



(45) 



where e| 3 is the permittivity component at constant strain 
under the assumption of plane stress. 

The normalized frequency response function (FRF) of the 
power output is 



P{co) = 



W\ 2 



(46) 



Rlg 2 Q) 2 

where g is the ground acceleration (9.81 m s" 2 ) and V is the 
amplitude of the steady state response in equations (40) and 
(41). 

Although the power output of a piezoelectric harvester 
depends on the circuit design, in this paper only the 
mechanical effect on the electromechanical behavior of an 
optimal bimorph piezoelectric system is considered. 

4. Numerical results 

In this section a case study is presented to show the 
performance of the device based on the analytical model given 
in section 3. The geometric and material properties of the 
bimorph cantilever investigated are shown in tables 1 and 2 
respectively. 



Table 1. Geometric properties of the bimorph cantilever. 
Piezoceramic Substructure 

Length (L c ) (mm) 100 100 

Width a (b) (mm) 15 15 

Thickness (h p , h $ ) (mm) 0.15 0.25 



a The width for a dual beam is 2 nm ; 
value that is used in the model. 



15 mm. This is the 



Table 2. Material properties of the bimorph cantilever. 





Piezoceramic 


Substructure 


Material 


PZT-5A 


Aluminum 


Elastic modulus (cf 1? F s ) (GPa) 


61 


70 


Mass density (p s , p p ) (kg m -3 ) 


7750 


2700 


Piezoelectric constant 


-10.4 


— 


(e 3 i)(Cm- 2 ) 






Permittivity constant 


13.3 


— 


(4)(nFm" 1 ) 







Table 3. Modal frequency of the beam without tip mass, stiffness or 
axial load. 



Mode 


12 3 4 5 


Frequency (Hz) 
Modal damping ratio f r 


29.91 187.47 524.96 1028.7 1700.54 
0.01 0.012 0.03 0.059 0.097 



Figure 6 displays the normalized displacements and 
strain mode shapes for the first three vibration modes of the 
cantilevered beam without tip stiffness or axial load. The 
natural frequencies for the first five modes are given in table 3. 
The FRF of the power output for a series connection can 
be obtained by equations (40) and (46). Figure 7 shows the 
FRFs of power outputs, normalized to the square of the base 
acceleration (i.e., g 2 ), for different load resistances. To select 
the optimal load resistance, the maximum of the power output 
for each load resistance was calculated and the results are 
shown in figure 8. It can be observed that there are two optimal 
power operating points. One of these maxima is related to 
the mechanical resonance whereas the other one is related 
to the unit resonance as a result of the piezoelectric circuit 
contributions to the overall stiffness and mechanical behavior 
of the device [34]. To investigate the behavior of the proposed 
device the second optimal point (i.e., 368.9 f^W g -2 , 109 k£2) 
is used throughout this paper. 

4.1. Parametric study 

Figure 9 shows the effects of the two parameters K and N, 
defined in section 3, on the natural frequency dimensionless 
parameter f n that is defined as the ratio between the frequency 
of the nth mode and the natural frequency of the first vibration 
mode without tip stiffness or axial load. Figure 9(a) shows the 
variation of the dimensionless natural frequencies of the first 
two vibration modes of the bimorph with tip stiffness (K), 
in the absence of axial load (N). As is expected, the natural 
frequencies of the first two modes increase as the stiffness 



Smart Mater. Struct. 21 (2012) 125025 



E Dehghan Niri and S Salamone 



a) 











/ 




Mode 1 






Mode 2 








Mode 3 




/ \ 














/ / 


\ 
\ 


"O^ 






\\ / 

'--,. .S "" 




\ 




Figure 6. Normalized (a) displacement and (b) strain mode shapes of the cantilevered beam without tip stiffness or axial load for the first 
vibration modes. 
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Figure 7. The power frequency response function of the bimorph without tip stiffness or axial load for series connection of the 
piezoceramic layers with (a) focus on the first three vibration modes and (b) focus on the first vibration mode. 
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Figure 8. Variation of the maximum of power output with load 
resistance for all excitation frequencies and considering the five 
modes of vibration for the cantilevered beam without tip stiffness or 
axial load and considering a series configuration for the electric 
circuit. 

of the dynamic system increases. It is noteworthy that the 
boundary conditions of the cantilevered harvester beam with 
a tip stiffness shift from clamped-free to clamped-pinned 



as K -¥ oo. In contrast, from figure 9(b) it can be inferred 
that, in absence of tip stiffness, by increasing the axial load 
the frequencies of the first two modes decrease. The reason 
behind this behavior is that an axial load provides a negative 
geometric stiffness, making the system softer. These two 
opposite effects can be combined by the proposed mechanism 
to control the dynamic response of the bimorph. 

Figure 10 shows the variation of the first dimensionless 
modal frequency parameter versus the dimensionless tip 
stiffness and axial load. It should be mentioned that the 
frequency range covered by the device can be modified by 
changing the device parameters, such as the spring stiffness 
K s . 

The figure shows that the resonance frequency can be 
tuned to a wide range of frequencies using a combination 
of tip axial load and stiffness. Although the ability to tune 
the harvester with tip stiffness and axial load to very low 
(close to zero in the border of the stable region) and very 
large frequencies (in the limit of the resonance frequency 
of the clamped-pinned beam) is proved, the performance of 
the harvester in power generation should be appropriate in 
comparison with the power generation of a harvester without 
tip stiffness or axial load. 

For evaluation of the efficiency of the bimorph with 
different tip stiffness and axial load combinations (P = 
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Figure 9. Variation of the dimensionless frequency parameter f n of the first and second vibration modes: (a) ratio of tip stiffness (no axial 
load), (b) ratio of tip axial load (no tip stiffness). 




Figure 10. 

of the first 
load. 



Variation of the dimensionless frequency parameter f n 
vibration mode for combination of tip stiffness and axial 




o o 



Figure 11. Variation of the dimensionless power generation 
parameter with respect to the dimensionless tip stiffness and axial 
load. 



1 when K = and N = 0), the dimensionless power 

aenerarirm narampfpr P — Pmax (with ti P stlffness and axial load ) 

generation parameter, r - Pmax (without tip stiffne ss and axial load) ' 
is shown in figure 1 1 . 

It can be seen that, in the absence of tip axial load, the 
power generation decreases by increasing the tip stiffness; 
this is caused by the increasing stiffness of the system and 
reduction in the displacement field of the harvester (reduction 
in the strain field as well). However, as can be observed 
in figure 11, the application of the axial load compensates 
the reduction of power generation even if the tip stiffness 
is increased. The physics behind this behavior is that the 
reduction in the general stiffness of the dynamic system as 
a result of the additional axial load leads to an increase of the 
displacement field of the harvester (the strain field increases). 
Hence, by applying the combination of tip stiffness and axial 
load, which is the main advantage of the proposed device, 
one can change the first modal frequency of the bimorph 
(figure 10) without significant effect on the power generation 
(figure 11). The specific design parameters of this device to 
tune to different frequencies and the corresponding numerical 
analysis are discussed in section 4.2. 



Table 4. The design properties of the device for the case study. 



(Note that ao = tan 



K 


N 





"max 


e (cm) 


a (cm) 


19.74 


1.051 


9.5° 


25.5° 


5 


2.5 



4.2. A design specification for passive tuning of the device to 
diffe rent frequencies 

Table 4 summarizes the design properties including the spring 
stiffness K s , spring initial preload N s o, initial spring angle 
$0, maximum spring angle £ max and e and ao shown in 
figure 2. These design specifications are used to investigate the 
performance of the device. However, it should be pointed out 
that these parameters must be designed to obtain the frequency 
range of interest. 

Figure 12 shows the relationship between d and 6, 
obtained from equations (1) and (2), for three different values 
of ao. Two main observations can be made: (1) the sensitivity 
of the angle to d can be modified by changing ao and (2) 
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Figure 12. Variation of the angle 9 versus the displacement d for 
different ao. 



small displacements (d) (e.g., between and 1.4 cm) can be 
used to cover a range of angle 9 between 9.5° and 25.5°. 

Figure 13 shows the power output FRF calculated in the 
range of angles covered using the design parameters in table 4. 
While the bimorph without tip mass or stiffness is just tuned to 
29.1 Hz, this mechanism allows an operator to manually tune 
the device in a range of frequencies between 6 and 61.87 Hz, 
just by changing d between and 1 .4 cm. 

Another interesting observation is that an initial increase 
of the maximum power output is followed by a decreasing 
trend. This behavior is caused by the strain distribution 
throughout the length of the beam, which affects the integrand 
in equation (44). Variations in the tip stiffness and axial 
load affect the natural frequencies and mode shapes; as a 
result, the modal strain along the length of the beam changes. 
Figure 14 shows the normalized first displacement and strain 
modal shape for different values of the angle 9. In the first 
stage, increasing the angle changes the modal strain in a way 
that increases the integrand (increasing trend in figure 13). 
After the power FRF has reached a certain value (or certain 



angle of the changing mechanism) this integrand decreases 
because of cancelation of the electrical output. A sign change 
for a vibration strain mode for a bimorph with a continuous 
electrode results in cancelation of the electrical output [33]. 
The change in strain sign for the first vibration mode is 
apparent in figure 14(b). Figure 15 shows that the effect of 
tip stiffness and axial load on the second strain mode is not as 
pronounced as on the first mode in this particular design. 

To evaluate the stability of the device, figure 16 shows 
the variation of the normalized equivalent tip stiffness and 
axial load in the range of angle 9 obtained using the design 
parameters in table 4. The compressive critical axial load for 
this case study without tip stiffness is TV = ^^ — 6.3471 N. 
The approximate stiffness of the cantilevered beam without 
axial load or tip stiffness, assuming a polynomial mode shape, 
is K — ^M — 77.17 N m _1 . It can be seen that change of 
the angle 9 within the range covered by the assumed design 
proprieties (e.g., 9.5°-25.5°) does not affect the stability of 
the device (i.e., the corresponding values of the equivalent tip 
stiffness and axial load stay within the stable region). It should 
be noted that one can use different design specifications listed 
in table 4 to change the combination of tip axial load and 
stiffness at changing mechanism d or angle 9 and thus match 
the desired design requirements path (the path is the inclined 
line shown in figure 16 for from 9 = 9.5° or d = to 9 = 
25.5° or d = 1.4 cm for the design specifications listed in 
table 4). 

The dimensionless power generation P of the system at 
different frequencies (different mechanism displacement d or 
angle 9) is shown in figure 17. As was discussed, although 
generally the power generation tends to decrease by increasing 
the tip stiffness, the tip axial load can compensate this effect 
reasonably. 

Figure 18 shows the dimensionless fundamental modal 
frequency of the device corresponding to different angles. 
This device with the stated design specifications can cover a 
frequency range between and almost two times the bimorph 
resonance frequency without tip stiffness and axial load. 

Indeed, the line location, inclination and length (from 
9 = 9.5° to 25.5°) shown in figures 16-18 can be changed 
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Figure 13. Power FRFs of the device for different values of the angle of the design specifications for (a) the three first modes and (b) the 
first mode. 
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Figure 14. Variation of the normalized (a) displacement and (b) strain mode shapes for the first vibration mode with the angle for the 
proposed device. 




Figure 15. Variation of the normalized (a) displacement and (b) strain mode shapes of the second vibration mode with the angle for the 
proposed device. 
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Figure 16. Equivalent tip stiffness and axial load at different angles 
for the design specifications. 

by selecting different design specifications like the spring 
stiffness K s , spring initial preload N s o and e to cover different 
ranges of frequencies with different paths. 
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Figure 17. Variation of the dimensionless power generation 
parameter with respect to the dimensionless tip stiffness and axial 
load of the device with the design specifications shown in table 4. 



To tune the frequency of the specified device to 
different frequencies one can adjust the displacement of 
the mechanism d according to figure 19. To provide the 
mechanism displacement d, it is suggested that the screw 
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Figure 18. Variation of the dimensionless frequency parameter f n of 
the first vibration mode for a combination of tip stiffness and axial 
load of the device with the design specifications shown in table 4. 




Figure 19. Variation of the tuning frequency with respect to the 
displacement of the mechanism d and the corresponding angle 9. 



connected to the simple joint at which the two levers meet 
each other should be used, figure 1 (this mechanical joint is 
horizontally and torsionally free). This screw can be found at 
the clamped connection at the end of the cantilevered energy 
harvester beam (figure 1). The pitch of this screw can be 



<* 



stable region 




Figure 21. The equivalent tip stiffness and axial load in the 
proximity of the unstable region for a few selected N CT of the design 
specifications. 



simply designed in such a way as to relate different screw 
revolutions to different frequencies. 

To evaluate the sensitivity of the resonance frequency 
of the device to the displacement (d) of the mechanism and 
spring angle two sensitivity criteria are used, 



s d - 



4A 
dd 



and 



So 



4/i 
d0' 



(47) 



These two criteria are the derivatives of the fundamental 
frequency (i.e., frequency of the first modal vibration) with the 
displacement of the mechanism, d, and angle, 9, respectively. 
Figure 20 shows variation of these two criteria with respect 
to the displacement of the mechanism, d, and angle, 0. For 
these specific design parameters the sensitivity of the system 
decreases as d and 6 increase. At the beginning id — or 6 = 
9.5°), the fundamental vibration frequency of the bimorph can 
be changed by 35 Hz by just a one millimeter displacement of 
the mechanism; this value reduces to 1 Hz at d equal to 1.4 cm 
(0 = 25.5°). 




Figure 20. The sensitivity of the tuning frequency f\ to (a) the displacement d of the mechanism and (b) the angle 9. 
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Frequency [Hz] Frequency [Hz] 

Figure 22. Power FRFs of the device for different N cr of the design specifications for (a) the two first modes and (b) the first mode. 




Figure 23. Variation of the normalized (a) displacement and (b) strain mode shapes of the first vibration mode for the proposed device with 
increasing N cr . 



4.3. Low frequency tuning 

One of the main advantages of this device, apart from its 
capability to be passively tuned to different frequencies with 
small adjustments of the sliding mechanism, is that it can 
be theoretically tuned to frequencies close to zero, as can be 
observed in figure 13. In fact, the flat FRF achieved for an 
angle of 6 =9.5° (i.e., frequencies smaller than 6 Hz) shows 
that this device does not require a large proof mass to be tuned 
to low frequencies. 

To show the capability of the device to provide flat and 
wideband power output FRFs for low frequencies, let us 
consider a set of normalized tip stiffnesses and axial loads 
in the borderline region, as shown in figure 21. Figure 22 
shows the corresponding power output FRFs. As the boundary 
condition of the bimorph with the axial load on the stable 
border gets closer to the clamped-pinned boundary condition, 
the maximum power output decreases and the device FRF 
supports a wider band of frequencies. As a result, the 
device can be tuned to a wider range of low frequencies in 
the borderline region, without losing significant efficiency. 
Further, the maximum power output for the second mode 
increases by moving up along this region. Figure 23 shows 
the variation of the normalized first displacement and strain 



modal shape for the corresponding points shown in figure 21. 
The integrand of the first strain mode decreases gradually and 
this is why the FRF of the power output decreases (the first 
peaks in figure 22). On the other hand, the power output of 
the second mode increases as a result of changing the second 
strain node, which results in less cancelation of power. The 
normalized displacement and strain of the second mode are 
displayed in figure 24. 

5. Conclusion 

A novel dual cantilevered piezoelectric energy harvester 
with variable tip stiffness and axial load was presented. 
A differential equation for an axially loaded cantilevered 
beam with tip stiffness, based on Euler-Bernoulli beam 
assumptions, was developed and the characteristic function 
for multimode vibration of the beam based on the boundary 
conditions was derived. The frequency response function of 
the power output was used to evaluate the performance of 
the system. The effects of dimensionless tip axial load and 
stiffness on the resonance frequency and power generation of 
the harvester have been discussed. Increase of the tip stiffness 
increases the resonance frequency; conversely, compressive 
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Figure 24. Variation of the normalized (a) displacement and (b) strain mode shapes of the second vibration mode for the proposed device 
with increasing N cr . 



axial load decreases the resonance frequency of the harvester. 
It was shown that, although increase of the tip stiffness causes 
a reduction in power generation as a result of reduction of the 
displacement field, this can be compensated by application of 
an axial load that reduces the general stiffness of the dynamic 
system. 

The behavior of the proposed device with design spec- 
ifications in the stable region was numerically investigated 
and it was shown that it is able to passively tune to different 
frequencies by small adjustments of the proposed sliding 
mechanism. Furthermore, the results showed that this device 
can be used to tune the energy harvester to very low 
frequencies without the need of a large proof mass. 
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